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4.  INTRODUCTION 


The  goal  of  this  project  is  to  construct  a  functioning  compact  positron  tomograph, 
whose  geometry  is  optimized  for  detecting  prostate  tumors  with  molecular  tracers  such  as 
[’’Cjcholine  (carbon-11  choline),  to  confirm  the  presence,  absence  or  progression  of 
prostate  cancer.  Compared  with  a  conventional  positron  emission  tomography  (PET) 
camera,  the  proposed  tomograph  has  (1)  the  same  spatial  resolution  in  the  region  of  the 
prostate;  (2)  improved  detection  efficiency;  and  (3)  significantly  lower  cost,  which  would 
increase  clinical  availability.  See  Figure  1  for  a  sketch  of  the  positron  tomograph, 
consisting  of  detector  banks,  patient  and  patient  bed.  See  Figure  2  for  photographs  of  a 
detector  module  and  partially-assembled  camera.  See  Appendix  1  for  3-D  drawings  of  the 
positron  tomograph  with  some  gantry  details.  See  Appendices  2  and  3  for  technical 
papers  presented  at  scientific  meetings. 

5.  BODY  OF  THE  PROGRESS  REPORT 

Below  are  the  list  of  tasks  as  proposed  in  our  original  Statement  of  Work.  A  summary 
of  our  progress  and  current  status  for  each  sub-task  follows. 

Task  1  Design  two  external  curved  detector  banks;  Months  1-10: 

a)  Determine  optimal  detector  bank  geometry  based  on  simulations  and  calculations. 

b)  Design  shielding. 

c)  Design  gantry. 

d)  Design  phantom  bed. 

Task  2  Construct  the  camera;  Months  3-26: 

a)  Construct  gantry. 

b)  Purchase  and  test  CTI ECAT  HR+  PET  detector  modules. 

c)  Configure  detector  modules  into  two  detector  banks. 

d)  Construct  and  test  event  processing,  coincidence  processing,  and  readout  electronics. 

e)  Configure  system  with  phantom  bed. 

Task  3  Develop  iterative  reconstruction  algorithms;  Months  2-31: 

a)  Develop  iterative  algorithms  using  Monte  Carlo  generated  data. 

b)  Optimize  algorithms  using  real  data. 

Task  4  Characterize  camera  performance  with  phantoms;  Months  27-36. 


Task  1(a):  Determine  optimal  detector  bank  geometry 

Based  on  our  computer  simulations  and  calculations,  we  have  optimized  and  finalized 
the  design  of  the  PET  camera  [1,2].  Figure  1  shows  the  transaxial  and  sagittal  views  of 
the  camera.  The  bottom  bank  is  fixed  below  the  patient  bed,  and  the  top  bank  moves 
upward  for  patient  access  and  downward  for  maximum  sensitivity.  Each  bank  consists  of 
two  axial  rows  of  20  CTI  ECAT  HR+  detector  modules,  forming  two  arcs  with  a  minor 
axis  of  45  cm  and  a  major  axis  of  70  cm.  Our  original  proposal  allowed  for  only  a  single 
axial  row,  but  patient  alignment  will  be  easier  with  greater  axial  coverage;  we  negotiated 
with  CPS  Innovations  to  receive  the  CTI  ECAT  HR-i-  detector  modules  at  below  cost  (to 
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afford  the  2nd  axial  row).  Individual  detector  modules  are  angled  to  point  towards  the 
camera  center  in  order  to  optimize  resolution  at  the  center  near  the  prostate  at  the  price  of 
some  resolution  degradation  at  the  edges.  Both  detector  banks  can  be  tilted  to  image  the 
prostate  while  minimizing  attenuation  (i.e.,  above  the  buttocks  and  below  the  stomach, 
see  Figure  1),  but  the  gantry  allows  zero  tilt  for  thin  patients. 

Task  1(b):  Design  shielding 

The  camera  design  includes  shielding  (see  Figure  lb).  Inter-plane  tungsten  septa 
extend  beyond  the  scintillator  crystals  to  reduce  background  events  from  random 
coincidences  and  from  photons  that  Compton  scatter  in  the  patient.  Lead  shields  are  also 
used  on  the  ends  to  reduce  activity  from  outside  the  field  of  view.  Our  original  septa 
design  consisted  of  septa  with  an  uniform  length  of  5  cm,  based  on  geometric 
calculations  with  reference  to  current  clinical  whole-body  scanners.  However,  we  have 
recently  used  Monte  Carlo  simulations  to  further  investigate  several  novel  septa  designs, 
in  addition  to  considering  manufacturing  techniques.  We  now  have  two  improved  septa 
designs  [3]  with  septa  that  vary  in  length  depending  on  the  arc  position.  The  first  design 
places  septa  (0.8  mm  thick,  with  a  length  that  varies  from  2  cm  at  the  minor  axis  to  5  cm 
at  the  major  axis)  axially  between  each  scintillator  crystal.  The  second  design  places 
septa  (2.4  mm  thick,  with  a  length  that  varies  from  10  cm  to  20  cm)  axially  between 
detector  blocks.  We’ll  choose  between  these  two  novel  designs  shortly,  based  on  cost  and 
expected  performance. 

Task  1(c);  Design  gantry 

The  gantry  and  electronics  housing  is  designed.  Appendix  1  shows  3-D  drawings  of 
the  camera  with  many  of  the  gantry  and  electronics  housing  details.  The  upper  bank  is 
mounted  on  a  main  aluminum  back-plate  in  order  to  control  the  vertical  motion  using  a 
hand-crank,  allowing  upward  movement  for  patient  access  and  downward  to  maximize 
sensitivity.  A  second  hand-crank  is  used  to  adjust  the  overall  tilt  of  both  upper  and  lower 
banks.  The  near  end  of  Appendix  1(a)  also  shows  the  electronics  housing;  see  Appendix 
2  for  further  electronics  details.  The  electronics  used  for  another  project  is  similar  to  this 
system  [4].  The  ±5V  and  high-voltage  power  supplies  are  visible  in  Appendix  lb. 

(a)  (b) 


Figure  1:  Positron  tomograph  for  prostate  imaging,  (a)  Drawing  of  atransaxial  view  through  prostate,  showing  the  patient 
centered  between  two  detector  banks  (shielding  not  shown).  The  individual  detector  modules  are  angled  to  point  towards  the 
prostate,  (b)  Drawing  of  the  sagittal  view  with  shielding  shown.  The  bottom  arc  is  fixed  below  the  patient  bed,  whereas  the 
top  arc  adjusts  vertically  for  patient  access  and  compactness.  Both  detector  banks  are  tilted  and  positioned  as  close  as 
possible  to  the  prostate,  which  improves  sensitivity  and  minimizes  attenuation. 
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Task  1(d):  Design  phantom  bed 

Our  proposal  is  limited  to  the  design,  construction  and  characterization  (using 
phantoms)  of  a  compact  prostate  imaging  system.  Since  human  studies  are  beyond  the 
scope  of  this  proposal,  we  have  only  designed  a  simple  phantom  bed  at  this  time. 


Task  2(a):  Construct  gantry 

The  gantry  and  electronics  housing  are  under  construction,  and  completion  is 
expected  before  April,  2004  {i.e.,  month  26).  All  gantry  components  have  been  machined 
or  purchased,  and  the  main  gantry  assembly  is  complete.  Assembly  of  the  electronics 
housing,  cabling,  and  power  supplies  is  in  progress.  Figure  2b  shows  a  photograph  of  the 
partially-assembled  camera,  taken  before  the  lead  shielding  was  glued  in  place. 

Task  2(b):  Purchase  and  test  CTI  ECAT  HR+  PET  detector  modules 

We  have  acquired  all  CTI  ECAT  HR+  detector  modules  and  electronics  boards  from 
CPS  Innovations  at  a  reduced  cost.  All  detector  modules  have  been  tested  and 
demonstrated  to  work  properly. 

Task  2(c):  Configure  detector  modules  into  two  detector  banks 

Figure  2  shows  photographs  of  a  detector  module  and  partially-assembled  camera.  All 
80  detector  modules  will  be  configured  into  two  detector  banks  as  soon  as  the  gantry  and 
electronics  housing  assembly  is  complete. 

Task  2(d):  Construct  and  test  event  processing,  coincidence  processing,  and  readout 
electronics 

We  have  tested  and  modified  the  CPS  Innovations  electronics  and  data  acquisition 
software  for  our  system,  as  well  as  designed  and  constructed  the  required  conversion 
boards  (for  the  Analog  Subsection  boards).  We  have  tested  the  event  processing, 
coincidence  processing  and  readout  electronics  using  test  pulses.  We  are  currently 
modi  fying  software  (e.g.,  software  to  calibrate  the  detector  modules). 

Task  2(e):  Configure  system  with  phantom  bed 


We  have  not  yet  begun  this  task.  We  expect  to  configure  the  system  with  the  phantom 
bed  this  year. 


Figure  2:  (a)  CTI  ECAT  HR+  (detector  moidule  with  an  aluminum  gui(ie  epoxietd  along  its  midline,  (b)  Photograph  of  the 
partially-assembled  camera  with  a  single  axial  row  of  detector  modules.  Each  detector  module  has  an  aluminum  guide  that 
slides  along  a  groove  in  the  arc  plate  and  is  stopped  (and  secured  against  the  arc  plate)  using  a  mortise  and  tenon.  The 
detector  modules  are  individually  angled  to  point  towards  the  center  of  the  camera. 
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Hgure  3:  (a)  Body  outline  of  patient  with  2.5  cm  diameter  circular  “prostate”  divided  into  quadrants  with  “tumors”  placed 
in  the  upper  right  and  lower  left  quadrants  (white  =  high  uptake),  assuming  a  tumor  to  background  ratio  of  2:1.  The 
“patient”  is  37.5  cm  wide  and  18.0  cm  deep,  (b)  Reconstructed  image  of  the  phantom  using  an  iterative  maximum 
likelihood  algorithm  that  assumes  a  30%  scatter  fraction  and  20%  randoms  fraction.  The  image  represents  745  kcounts  in  a 
single  slice,  which  should  be  achievable  with  a  6-minute  scan  after  a  10  mCi  injection.  Both  tumors  are  clearly  visible. 


Task  3(a):  Develop  iterative  algorithms  using  Monte  Carlo  generated  data 

The  reconstruction  algorithm  is  under  development.  We  have  done  a  significant 
amount  of  simulation  work  to  evaluate  the  expected  camera  performance  as  well  as  to 
optimize  the  camera  design.  We  have  determined  the  size  of  the  region  with  high  spatial 
resolution  within  our  PET  camera,  showing  that  there  will  be  good  resolution  in  the  10 
cm  diameter  central  region  near  the  prostate  for  a  wide  range  of  patient  sizes  [2],  We 
have  addressed  the  issue  of  irregular  and  incomplete  sampling  due  to  the  side  gaps 
between  the  detector  banks  [1]. 

Monte  Carlo  methods  have  been  used  to  study  the  reconstruction  of  extended  prostate 
sources.  Reconstructions  were  performed  with  an  iterative  maximum  likelihood 
technique  that  assumes  statistical  noise  representing  745  kcounts/slice,  a  tumor  activity 
concentration  that  is  two  times  the  background  activity  concentration,  30%  scatter 
fraction,  and  20%  randoms  fraction.  The  reconstructions  were  then  post-filtered  with  a 
very  conservative  Gaussian  function  of  sigma  =  1  pixel  =  2  mm.  We  have  demonstrated 
that  we  can  reconstruct  nearly  artifact-free  images,  as  well  as  detect  and  differentiate 
partial  and  whole  prostate  tumors,  using  Monte  Carlo  data  equivalent  to  a  6-minute  scan 
after  a  10  mCi  radiopharmaceutical  injection  (see  Figure  3).  The  expected  camera 
performance  is  discussed  in  further  detail  in  Appendix  2.  We  have  explored  the  necessary 
data  corrections,  such  as  attenuation  correction  and  normalization.  Further  development 
is  still  needed,  including  the  incoiporation  of  scatter  and  the  optimization  of  the  algorithm 
using  Monte  Carlo  data. 

Task  3(b):  Optimize  algorithms  using  real  data 

We  have  not  yet  begun  this  task.  We  expect  to  optimize  the  algorithm  using  real  data 
this  year. 

Task  4:  Characterize  camera  performance  with  phantoms 

We  have  not  yet  begun  this  task.  Characterization  of  camera  performance  with 
phantoms  will  be  described  in  the  next  progress  report. 
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6.  KEY  RESEARCH  ACCOMPLISHMENTS 


•  The  camera  design  was  optimized  and  completed. 

•  All  gantry  components  were  fabricated,  and  the  main  gantry  assembly  was  completed. 

•  Designs  for  electronics,  packaging,  mounting  and  interconnects  were  completed. 
Construction  will  be  completed  before  April,  2004. 

•  All  detector  modules  and  electronics  were  acquired  and  tested. 

•  Extensive  simulations  were  used  to  evaluate  the  expected  camera  performance. 

7.  REPORTABLE  OUTCOMES 

J.  S.  Huber,  J.  Qi,  S.  E.  Derenzo,  et  al.,  “The  Development  of  a  Compact  Positron 
Tomograph  for  Prostate  Imaging,”  LBNL-50790,  2002.  Oral  presentation  at  the  2002 
IEEE  Medical  Imaging  Conference. 

J.  Qi,  J.  S.  Huber,  R.H.  Huesman,  W.  W.  Moses,  S.  E.  Derenzo,  and  T.  F.  Budinger, 
“Septa  Design  for  a  Prostate-Specific  PET  Camera,”  IEEE  Trans  Nucl  Sci,  Submitted, 
2003.  Oral  presentation  at  the  2003  IEEE  Medical  Imaging  Conference. 

8.  CONCLUSIONS 

Task  1  has  been  completed.  Design  of  the  detector  bank  geometry,  electronics, 
gantry,  electronics  housing  and  patient  bed  was  optimized  and  completed.  Septa  design 
was  improved. 

Task  2  is  in  progress  and  on  schedule.  The  gantry  components  were  fabricated,  and 
the  main  gantry  assembly  was  completed.  All  detector  modules  and  electronics  were 
acquired  and  tested.  DAQ  software  was  modified,  and  further  software  development  is  in 
progress.  A  phantom  bed  was  constructed.  The  system  will  be  configured  once  the  gantry 
and  electronics  housing  assembly  is  complete.  No  problems  related  to  this  task  are 
currently  anticipated. 

Task  3  is  in  progress  and  on  schedule.  Extensive  simulations  were  used  to  evaluate 
the  expected  camera  performance,  demonstrating  that  our  PET  camera  design  will  image 
prostate  tumors  with  good  spatial  resolution  and  image  contrast  at  low  cost  relative  to 
commercial  PET  scanners,  in  agreement  with  our  original  proposal.  A  preliminary 
iterative  reconstruction  algorithm  was  developed,  but  further  development  is  necessary. 
Task  4  has  not  yet  begun. 
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10.  APPENDICES 


Appendix  1:  3-D  color  drawings  of  the  PET  camera  with  some  gantry  details;  (a)  as 
viewed  from  the  electronics  end  and  (b)  as  viewed  from  the  power  supply/doctor  access 
end. 

Appendix  2:  J.  S.  Huber,  J.  Qi,  S.  E.  Derenzo,  et  al.,  “The  Development  of  a  Compact 
Positron  Tomograph  for  Prostate  Imaging,”  LBNL-50790,  2002. 

Appendix  3:  Preprint  of  J.  Qi,  J.  S.  Huber,  R.  H.  Huesman,  et  al.,  “Septa  Design  for  a 
Prostate-Specific  feT  Camera,”  IEEE  Trans  Nucl  Sci,  Submitted,  2003. 
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Back-jalate 


Lead  Shielding 


Appendix  1:  3-D  color  drawings  of  the  prostate  PET  camera  with  some  gantry  details.  Only  a  few 
detector  modules  are  shown  in  place.  The  arc  plates  are  mounted  on  thicker  aluminum  arc  back-plates  that 
also  hold  the  lead  shielding  in  place.  The  arc  plates  (pink),  arc  back-plates  (blue)  and  lead  shielding  (green) 
are  symmetrical  since  there  are  two  axial  rows  of  detector  modules.  The  upper  bank  is  mounted  on  a  main 
aluminum  back-plate  in  order  to  control  the  vertical  motion  using  a  hand-crank,  allowing  upward 
movement  for  patient  access  and  downward  to  maximize  sensitivity.  A  second  hand-crank  is  used  to  adjust 
the  overall  tilt  of  both  upper  and  lower  banks.  Most  of  the  cabling  is  not  shown,  as  well  as  the  patient  bed. 
(a)  The  near  end  of  the  drawing  shows  the  electronics.  The  six  Detector  Head  Interface  boards  are  mounted 
horizontally  to  the  housings  that  hold  their  corresponding  Analog  Subsection  boards.  The  Coincidence 
Processor  board  is  mounted  on  the  cabinet  door,  (b)  The  near  end  of  the  drawing  shows  the  power  supplies. 
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The  Development  of  a  Compact  Positron 
Tomograph  for  Prostate  Imaging 

J.  S.  Huber,  Member,  IEEE,  J.  Qi,  Member,  IEEE,  S.  E.  Derenzo,  Fellow,  IEEE,  W.  W.  Moses, 
Senior  Member,  IEEE,  R.  H.  Huesman,  Fellow,  IEEE,  and  T.  F.  Budinger,  Member,  IEEE 


Abstract — ^We  give  design  details  and  expected  image 
results  of  a  compact  positron  tomograph  designed  for  prostate 
imaging  that  centers  a  patient  between  a  pair  of  external 
curved  detector  banks  (ellipse:  45  cm  minor,  70  cm  major 
axis).  The  bottom  bank  is  fixed  below  the  patient  bed,  and  the 
top  bank  moves  upward  for  patient  access  and  downward  for 
maximum  sensitivity.  Each  bank  is  composed  of  two  rows 
(axially)  of  20  CTI  PET  Systems  HR+  block  detectors,  forming 
two  arcs  that  can  be  tilted  to  minimize  attenuation.  Compared 
to  a  conventional  PET  system,  our  camera  uses  about  one- 
quarter  the  number  of  detectors  and  has  almost  two  times 
higher  solid  angle  coverage  for  a  central  point  source,  because 
the  detectors  are  close  to  the  patient.  The  detectors  are  read  out 
by  modified  CTI  HRRT  data  acquisition  electronics.  The 
individual  detectors  are  angled  in  the  plane  to  point  towards 
the  prostate  to  minimize  resolution  degradation  in  that  region. 
Inter-plane  septa  extend  5  cm  beyond  the  scintillator  crystals 
to  reduce  random  and  scatter  backgrounds.  Average-  to  large- 
size  patients  will  not  be  fully  encircled  by  detector  rings, 
causing  incomplete  sampling  due  to  the  side  gaps.  Monte  Carlo 
simulations  (including  20%  randoms  and  30%  scatter 
fractions)  demonstrate  the  feasibility  of  detecting  and 
differentiating  prostate  tumors  with  a  tumor  to  background 
ratio  of  2:1,  using  a  number  of  counts  that  should  be 
achievable  with  a  6  minute  scan  after  a  10  mCi  injection  {e.g., 
carbon-11  choline).  Simulations  also  predict  minimal  blurring 
in  the  10  cm  diameter  central  region  for  a  wide  range  of  patient 
sizes. 

I.  Introduction 

E  present  the  development  of  a  compact  positron 
tomograph  optimized  to  image  the  prostate.  This 
instrument  images  radiopharmaceuticals  that  specifically 
localize  in  the  prostate  to  confirm  the  presence,  absence  or 
progression  of  disease.  It  has  approximately  four  times  fewer 
detectors  than  a  conventional  whole-body  positron  emission 
tomograph  (PET),  which  will  reduce  the  cost  and  increase 
clinical  availability. 

Prostate  cancer  has  a  prevalence  and  diagnostic  rate  similar 
to  breast  cancer,  with  360,000  new  cases  diagnosed  each  year 
and  two  million  men  affected  by  the  disease  in  the  United 
States.  Prostate  cancer  suspicion  is  typically  based  on  an 
elevated  prostate-specific  antigen  (PSA)  level  or  a  suspicious 
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node  found  during  a  digital  rectal  exam.  Serum  PSA  values 
do  not  always  correlate  well  with  clinical  diagnosis  or 
outcomes  [1-3].  Palpation  is  subjective,  insensitive  and 
inexact;  more  than  half  of  all  cancers  detected  today  are  not 
palpable.  Treatment  decision  is  based  mainly  on  biopsy 
confirmation  of  prostate  cancer,  but  the  diagnostic  accuracy  of 
prostate  biopsies  is  problematic.  Typical  treatment  options 
include  radical  prostatectomy,  external  beam  irradiation, 
brachytherapy  (interstitial  implantation  of  radioactive  seeds), 
androgen  ablation  (hormone)  therapy,  or  "watchful  waiting." 
A  major  problem  with  prostatic  cancer  therapy  is  the  question 
of  when  to  treat  or  whether  to  treat  at  all.  This  is  particularly 
problematic  in  the  case  of  an  increased  PSA  level  with  non¬ 
diagnostic  repeated  biopsies  or  after  a  prostatectomy.  A  new 
imaging  technology  for  sensitive  detection  of  early  stage 
prostate  cancer  is  needed  to  confirm  initial  diagnosis  and  help 
guide  treatment  decisions.  In  addition,  a  new  method  is 
needed  to  assess  response  shortly  after  treatment  intervention. 

In  order  to  help  meet  these  needs,  we  are  building  a 
compact  PET  camera  optimized  to  image  the  prostate. 
Functional  PET  imaging  will  help  detect  malignant  tumors 
in  the  prostate  and/or  prostate  bed  to  confirm  an  elevated 
PSA  level,  as  well  as  possibly  help  determine  tumor 
"aggressiveness"  based  on  metabolic  uptake  levels,  in  order 
to  help  guide  whether  to  treat  suspected  prostatic  cancer. 
Although  not  optimized  to  detect  distant  metastatic  disease, 
this  compact  PET  camera  should  also  image  local  spread 
beyond  the  prostate  bed  to  help  guide  treatment  decisions 
such  as  whether  a  narrow  or  wide  irradiation  treatment  field  is 
needed. 

Promising  PET  radiopharmaceuticals  have  recently 
demonstrated  outstanding  results  in  the  sensitive  detection  of 
prostate  cancer,  inspiring  a  new  interest  in  using  PET  for 
prostate  cancer  imaging.  Consistent  with  the  evidence  of 
increased  pool  size  of  choline  in  prostate  cancer  [4,  5],  Kara 
and  co-workers  demonstrated  prostate  tumor  uptake  by 
[*^C]choline.  They  find  that:  [”C]choline  clears  the  blood 
faster  than  FDG;  its  uptake  in  prostate  tumors  is  significantly 
higher  than  in  normal  and  surrounding  tissues  [6,  7], 
providing  excellent  tumor/normal  contrast;  and  bladder 
accumulation  is  minimal  if  the  correct  time  course  is  chosen 
[8]  which  is  a  major  advantage  over  FDG.  Additionally,  they 
observe  that  ["Cjcholine  PET  is  more  sensitive  for  detecting 
bone  metastases  than  bone  scintigraphy.  Therefore, 
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[^^CJcholine  is  an  attractive  PET  tracer  for  imaging  primary 
and  metastatic  tumors  of  the  prostate  and  potentially  for  other 
regions  of  the  body  [9-13].  Fig.  1  shows  a  [*^C]choline 
image  of  prostate  cancer  before  and  after  therapy, 
demonstrating  the  ability  to  detect  prostate  carcinoma  using 
choline  and  possibly  its  analogs.  Other  "C 

radiopharmaceuticals  are  also  under  investigation  for  prostate 
cancer  PET  imaging,  including  ["Cjacetate  and 
[^^Cjmethionine. 

imaging  has  the  advantage  of  a  longer  half-life,  which 
increases  clinical  viability  since  an  on-site  cyclotron  facility 
would  not  be  necessary  (as  it  would  for  ’’C  imaging).  Many 
groups  have  shown  that  PET  imaging  with 

[**F]fluorodeoxyglucose  (FDG)  is  not  a  good  technique  for 
prostate  cancer  diagnosis,  because  FDG  is  not  very  prostate 
specific  (with  a  SUV  of  typically  ~2)  and  the  uptake  is 
overcome  by  background  primarily  due  to  bladder 
accumulation.  The  SUV  (standard  uptake  value)  is  the  ratio 
of  activity  in  the  target  tissue  to  the  average  activity  in  the 
body.  However,  there  are  several  other  ‘®F 
radiopharmaceuticals  currently  under  investigation  for 
prostate  cancer  imaging,  including  [^^Fjfluorocholine  (FCH) 
[14,  15].  PET  images  using  [*^F]fluorocholine  demonstrate 
high  standardized  uptake  values  (eg.,  SUV  of  8),  indicating 
that  FCH  is  well  localized  in  the  prostate  cancer  and  can  be 
imaged  with  good  resolution  if  a  short  scan  time  (--five 
minutes)  is  used.  Therefore,  [^®F]fluorocholine  is  also  an 
attractive  PET  tracer  for  imaging  primary  and  metastatic 
tumors  of  the  prostate. 

II.  Camera  Design 

A.  Overview 

These  new  prostate  tracers  have  motivated  us  to  build  a 
low  cost  PET  camera  optimized  to  image  the  prostate. 
Coincidence  imaging  of  positron  emitters  is  achieved  using  a 
pair  of  external  curved  detector  banks,  one  placed  above  and 
one  below  the  patient.  Fig.  2  shows  the  transaxial  and 
sagittal  views  of  the  camera.  The  bottom  bank  is  fixed  below 
the  patient  bed,  and  the  top  bank  moves  upward  for  patient 
access  and  downward  for  maximum  sensitivity.  Each  bank 
consists  of  two  axial  rows  of  20  CTI  ECAT  HR+  block 
detectors,  forming  two  arcs  with  a  minor  axis  of  45  cm  and 
major  axis  of  70  cm.  Our  prostate  camera  has  about  one- 
fourth  the  number  of  detectors  as  in  a  conventional  PET 
system  because:  (a)  the  patient  is  not  fully  encircled  in  2D, 
(b)  an  elliptical  shape  is  used  instead  of  a  circular  one,  and  (c) 
the  axial  coverage  is  only  8  cm.  However,  since  the  average 
distance  to  the  detectors  is  approximately  one-half  that  of  a 
conventional  2D  PET  system,  the  solid  angle  for  a  central 
source  is  approximately  double  for  average-sized  men.  Thus, 
we  expect  to  achieve  improved  detection  efficiency  at  a  lower 
cost. 

Individual  detector  modules  are  angled  to  point  towards 
the  camera  center  near  the  prostate  location  in  order  to  reduce 
penetration  effects  (in  the  detector)  for  annihilation  photons 
originating  in  the  prostate.  This  increases  resolution 


selectively  in  the  region  of  the  prostate  and  is  an  unique 
feature  for  a  non-circular  camera  geometry.  Annihilation 
photons  fi*om  other  parts  of  the  field  of  view  (FOV)  will 
suffer  increased  penetration  effects,  but  these  FOV  regions  are 
less  important.  Both  detector  banks  can  be  tilted  to  image  the 
prostate  while  minimizing  attenuation  (i.e.,  above  the 
buttocks  and  below  the  stomach,  see  Fig.  2),  but  the  gantry 
allows  zero  tilt  for  thin  patients.  A  patient  of  average  size  is 
not  folly  encircled  in  2D,  which  results  in  irregular  and 
incomplete  sampling  due  to  the  side  gaps.  Despite  this 
incomplete  sampling,  we  are  able  to  reconstruct  nearly  artifact 
fi^ee  images  in  the  region  of  interest  by  using  an  iterative 
reconstruction  algorithm  [16]. 

The  camera  design  also  includes  shielding  [16].  Inter-plane 
septa  extend  5  cm  beyond  the  scintillator  crystals  to  reduce 
background  events  from  random  coincidences  and  fi-om 
photons  that  Compton  scatter  in  the  patient.  Lead  shields  are 
also  used  on  the  ends  to  reduce  activity  fi-om  outside  the  field 
of  view. 

B.  Electronics 

Our  camera  uses  modified  commercial  components  in  a 
novel  geometry.  We  are  using  80  CTI  ECAT  HR+  block 
detectors  that  are  three  attenuation  lengths  thick  for  good 
detection  efficiency  with  narrow  detector  elements  (/.e.,  8x8 
arrays  of  4.5  x  4.5  x  30  mm^  BGO  crystals)  to  achieve  good 
spatial  resolution.  We  are  using  modified  fi-ont  end, 
coincidence,  and  readout  electronics  developed  by  CTI  for  the 
HRRT  brain  imaging  PET  camera,  as  shown  in  Fig.  3.  Since 
we  are  creating  a  non-standard  camera  with  HR+  and  HRRT 
components,  a  custom  conversion  board  is  necessary.  Events 
are  detected  and  assigned  an  arrival  time,  the  crystal  of 
interaction  is  identified,  energy  qualification  is  performed, 
and  a  digital  word  is  formed  using  28  CTI  HRRT  Analog 
Subsection  boards.  The  output  signals  fi-om  the  Analog 
Subsection  boards  are  then  multiplexed  by  six  custom 
Detector  Head  Interface  (DHI)  boards,  which  are  based  on  the 
CTI  HRRT  DHI  design  but  each  services  a  maximum  of  1 5 
(rather  than  117)  detector  modules.  We  are  using  six  DHI 
boards  to  allow  coincidences  between  detector  modules 
within  the  same  bank  and  plane.  A  CTI  Coincidence 
Processor  identifies  singles  events  in  the  different  DHI  boards 
that  are  in  coincidence,  and  the  output  is  sent  to  a  Pentium- 
based  computer  system  that  accumulates  the  coincident  data, 
reconstructs  the  images,  and  displays  them. 

We  previously  reported  on  the  conceptual  design  of  this 
prostate  camera  [16].  We  have  moved  from  the  idea  stage  into 
the  building  stage.  We  have  acquired  all  HR+  detector  blocks 
and  electronics  boards.  We  are  currently  testing  the  CTI 
electronics  and  modifying  the  DAQ  software  for  our  system. 

C.  Gantry 

The  gantry  and  electronic  housing  is  designed  and  imder 
construction.  We  still  need  to' design  and  build  the  septa  and 
patient  bed.  Fig.  4  shows  a  photograph  of  the  top  and 
bottom  aluminum  arc  plates  with  some  of  the  CTI  ECAT  - 
HR+  detector  modules  mounted  in  place.  The  detector 
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modules  are  individually  positioned  to  point  towards  the 
center  of  the  camera,  as  discussed  in  Section  IIA. 

Fig.  5  shows  a  3-D  drawing  of  the  camera  with  many  of 
the  gantry  and  electronic  housing  details.  The  upper  bank  is 
mounted  onto  a  main  aluminum  back-plate  in  order  to 
control  the  vertical  motion  using  a  hand-crank,  allowing 
upward  movement  for  patient  access  and  downward  to 
maximize  sensitivity.  A  second  hand-crank  is  used  to  adjust 
the  overall  tilt  of  both  upper  and  lower  banks.  The  near  end 
of  Fig.  5  also  shows  the  electronic  housing.  The  six  DHI 
boards  are  mounted  horizontally  to  the  housing  that  holds 
their  corresponding  Analog  Subsection  boards.  The 
Coincidence  Processor  board  is  mounted  onto  the  cabinet 
door  for  compactness.  The  ±5V  and  high  voltage  power 
supplies  are  mounted  on  the  far  end. 

III.  Expected  Performance 
A.  Point  Source  Grid 

To  determine  the  size  of  the  region  of  high  spatial 
resolution,  we  have  performed  statistical  noise-free 
simulations.  Fig.  6  shows  reconstructed  point  sources  in  a 
grid.  Data  are  generated  with  interactions  throughout  the 
crystal,  assuming  an  attenuation  length  of  10  mm  in  the 
BGO  crystal.  A  preliminary  iterative  maximum  likelihood 
algorithm  is  used  in  the  reconstruction  without  correcting  for 
detector  penetration.  The  point  sources  are  5  cm  apart,  and 
the  grid  covers  the  majority  of  the  field  of  view.  All  point 
sources  are  clearly  visible  in  the  image  (black  =  high  uptake). 
As  expected,  increased  radial  blurring  due  to  penetration 
effects  is  seen  in  the  region  representing  the  outer  edge  of  the 
patient.  Since  we  are  optimizing  resolution  at  the  center  near 
the  prostate  at  the  price  of  resolution  at  the  edges,  distant 
metastases  will  suffer  from  some  resolution  degradation. 
However,  good  resolution  is  observed  in  the  10  cm  diameter 
central  region  near  the  prostate.  In  Fig.  6,  the  central  point 
has  both  a  radial  and  tangential  fwhm  of  2  mm.  More 
representative  of  the  expected  spatial  resolution,  the 
surrounding  eight  points  have  a  radial  fwhm  of  4  mm  and 
tangential  fwhm  of  3  mm.  Improved  resolution  for  more 
distant  points  is  expected  when  the  reconstruction  corrects  for 
detector  penetration,  although  this  process  increases  the 
statistical  noise  in  the  reconstruction. 

Assuming  that  the  detector  arcs  are  rigid  with  a  fixed  focal 
spot  for  each  arc,  both  focal  spots  are  at  the  center  for  the 
‘nominal’  gap  separation.  The  above  simulation  (Fig.  6) 
assumed  this  ‘nominal’  gap  distance.  In  order  to  investigate 
the  effect  of  varying  the  gap  between  the  arcs,  the  simulation 
is  repeated  for  a  gap  range  of  ±5  cm  when  the  focal  spots  are 
no  longer  in  the  center.  The  same  minimal  blurring  within  a 
10  cm  diameter  central  region  is  observed  for  both  cases, 
indicating  that  a  wide  range  in  patient  size  will  not 
compromise  the  resolution  near  the  prostate.  Fig.  7  shows 
the  reconstructed  point  source  grid  with  the  gap  distance 
increased  by  +5  cm  from  ‘nominal’  position. 

The  above  results  show  that  we  can  achieve  similar 
resolution  in  the  10  cm  diameter  center  region  as  an  ECAT 


HR+  scanner.  Using  the  nominal  configuration,  we  also 
compute  the  signal  to  noise  ratio  (SNR)  of  an  ideal  observer 
for  detecting  a  small  lesion  at  the  grid  locations.  The 
background  is  modeled  as  an  imiform  elliptical  region  (major 
axis  =  40  cm,  minor  =  30  cm).  For  the  central  point,  the 
resulting  SNR  of  the  new  system  is  about  27%  higher  than 
that  of  an  ECAT  HR+  scanner  under  the  same  situation,  and 
the  SNR  is  about  20%  higher  for  the  surrounding  eight 
points.  This  indicates  that  while  the  resolution  is  similar,  the 
increase  in  sensitivity  can  still  improve  lesion  detection.  The 
above  computation  is  for  one  2D  direct  plane.  No  random  or 
scattered  events  are  modeled  for  either  scanner. 

B.  Extended  Sources 

Monte  Carlo  methods  are  used  to  study  the  reconstruction 
of  extended  sources.  Reconstructions  are  performed  with  an 
iterative  maximum  likelihood  technique  that  assumes 
statistical  noise  representing  745  kcounts/slice,  a  tumor 
activity  concentration  that  is  two  times  the  background 
activity  concentration,  30%  scatter  fraction,  and  20% 
randoms  fraction.  The  reconstructions  are  then  post-filtered 
with  a  very  conservative  Gaussian  function  of  sigma  =  1 
pixel  =  2  mm.  We  simulate  self-attenuation  by  assuming  a 
constant  attenuation  coefficient  of  0.0095/mm  inside  the 
body.  Coincidences  between  detector  modules  within  the 
same  bank  and  plane  are  allowed.  Fig.  8(a)  shows  the  body 
outline  of  a  patient  with  an  extended  “prostate”  source  (2.5 
cm  diameter  circle)  with  “tumors”  in  the  upper  right  and 
lower  left  quadrants.  Fig.  8(b)  shows  a  reconstructed  image 
of  the  phantom  with  “tumors”  in  the  upper  right  and  lower 
left  quadrants  clearly  visible  (white  =  high  uptake), 
demonstrating  that  we  can  reconstruct  nearly  artifact-free 
images  as  well  as  detect  and  differentiate  partial  and  whole 
prostate  tumors.  Our  preliminary  rate  calculations  indicate 
that  745  kcounts  per  imaging  plane  will  be  achievable  with  a 
six  minute  scan  after  a  10  mCi  injection  [16].  A  relatively 
short  value  of  six  minutes  is  chosen  to  demonstrate  that  good 
images  can  be  obtained  even  if  the  bladder  fills  rapidly.  All 
published  ["CJcholine  or  [‘^FJfluorocholine  studies  done  to 
date  have  used  three  minutes  or  longer  (up  to  30  minutes) 
imaging  times. 

C  Data  Correction 

Our  design  has  sufficient  room  to  incorporate  a  single 
photon  or  positron  emitting  source  for  transmission  scans, 
but  it  would  add  to  the  camera  cost.  We  plan  to  correct 
without  an  additional  source  using  calculated  attenuation 
coefficients  based  on  body  contours  and  an  uniform 
attenuation  coefficient  [16].  Anatomical  boundaries  can  be 
obtained  from  the  outer  edges  of  emission  sinograms 
acquired  from  transverse  sections  [17].  It  may  also  be 
possible  to  use  heterogeneous  attenuation  coefficients  by 
identifying  tissue  types  (e.g.,  soft  tissue,  bone,  and  air)  based 
on  x-ray  CT  images. 

Because  of  the  unusual  geometry  of  our  camera,  the 
sensitivity  of  each  line  of  response  depends  on  the  distance 
between  the  two  detectors  and  the  face  angles  of  the  two 
detectors;  the  shorter  the  distance  is,  the  higher  the 
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sensitivity.  Normalization  will  include  geometric  correction 
(sensitivity  changes  caused  by  the  septa  that  are  not  modeled 
in  the  solid  angle  computation),  block  effect  correction, 
individual  detector  sensitivity  correction,  and  dead  time 
correction.  These  factors  will  be  obtained  using  a  scan  of  an 
uniform  activity  plane  source  such  as  used  in  conventional 
PET  scanners.  However,  because  of  our  specialized  geometry 
and  lack  of  rotational  symmetry,  the  method  for  deriving  die 
first  two  factors  will  be  different  from  that  of  circular  systems 
and  may  require  a  longer  scan  time  to  achieve  good  count 
statistics. 

Random  events  will  be  estimated  using  the  delayed 
window  technique.  Instead  of  pre-subtracting  histogrammed 
data,  delayed  events  will  be  saved  in  the  list  mode  data 
stream  with  a  tag  to  distinguish  them  from  prompt  events. 
For  our  statistical  reconstruction  algorithm,  the  estimated 
randoms  will  be  included  in  the  forward  model  in 
reconstruction. 

IV.  Discussion 

This  project  is  in  its  initial  building  phase,  so  there  are 
many  remaining  uncertainties,  including  which 
radiopharmaceutical  will  be  used,  the  specific  requirements 
for  the  radiopharmaceutical  and  imaging  procedure,  and  the 
accuracy  that  physicians  would  require  to  perform  effective 
diagnosis  and  evaluate  response  to  therapy.  However,  new 
PET  radiopharmaceuticals  have  recently  demonstrated 
promising  results  in  the  sensitive  detection  of  prostate  cancer. 
This  PET  camera  design  is  also  likely  to  out-perform  whole 
body  PET  for  this  task,  as  the  higher  sensitivity  will  help  for 
the  short  scans  that  appear  to  be  necessary  with  choline-like 
radiopharmaceuticals. 

Accurate  patient  positioning  is  also  critical  due  to  the 
limited  (8  cm)  axial  extent  of  the  camera.  A  manually  placed 
external  or  trans-rectal  positron-emitting  source  could  be  used 
to  position  the  prostate  near  the  center  of  the  FOV  of  the 
camera.  Alternatively,  the  prostate  is  visible  in  ultrasound 
images  acquired  with  an  external  transducer,  and  these  images 
could  be  used  to  position  the  patient  bed. 

Finally,  as  Fig.  1  shows,  it  is  difficult  to  identify  features 
other  than  the  prostate  (and  possibly  the  bladder)  in  the  PET 
images.  Local  anatomic  information  is  highly  desirable  to 
determine  the  location  of  the  disease  within  the  prostate,  as 
well  as  to  determine  whether  the  disease  has  spread  to  the 
prostate  bed.  This  information  could  be  obtained  with  dual 
modality  imaging.  X-ray  CT  provides  excellent  images  of  the 
abdomen,  but  the  cost  is  relatively  high.  Ultrasound  imaging 
with  a  trans-rectal  probe  provides  reasonable  detail  in  the 
region  of  the  prostate  and  the  cost  of  an  ultrasound  unit  is 
significantly  less  than  that  of  a  CT  imager,  so  we  plan  to 
explore  adding  co-registered  ultrasound  imaging  capability. 

V.  Conclusion 

Promising  new  PET  tracers  for  prostate  cancer,  such  as 
[^^Cjcholine  and  [^*F]fluorocholine,  have  motivated  us  to 
design  and  build  a  low  cost  PET  camera  optimized  for 
prostate  imaging.  Monte  Carlo  simulations  demonstrate  the 


feasibility  of  detecting  and  differentiating  partial  and  whole 
prostate  tumors  with  a  tumor  to  background  ratio  of  2:1, 
using  a  six  minute  scan  after  a  10  mCi  injection. 
Simulations  also  predict  minimal  blurring  in  the  10  cm 
diameter  central  region  for  a  wide  range  in  patient  size.  Thus, 
we  are  currently  building  a  camera  that  will  image  prostate 
tumors  with  good  spatial  resolution  and  image  contrast  at 
low  cost  relative  to  commercial  PET  scanners. 
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Fig.  3:  Diagram  of  the  complete  data  collection  chain,  showing  the  custom 
LBNL  and  CTI  components  and  their  inter-relations.  The  camera  will  use 
80  CTI  ECAT  HR+  block  detectors,  28  LBNL  conversion  boards,  28  CTI 
HRRT  Analog  Subsection  boards,  6  LBNL  custom  Detector  Head 
Interface  boards,  1  CTI  Coincidence  Processor  and  a  PC. 


Fig.  2;  Positron  tomograph  for  prostate  imaging,  (a)  Drawing  of  a  transaxial 
view  through  prostate,  showing  the  patient  centered  between  two  detector 
banks.  The  individual  detector  modules  are  angled  to  point  towards  the 
prostate,  (b)  Drawing  of  the  sagittal  view.  The  bottom  arc  is  fixed  below 
the  patient  bed,  whereas  the  top  arc  adjusts  vertically  for  patient  access  and 
compactness.  Both  detector  banks  are  tilted  and  positioned  as  close  as 


Fig.  4:  (a)  CTI  ECAT  HR+  detector  module  with  an  aluminum  guide 
epoxied  along  its  midline,  (b)  Photograph  of  the  top  and  bottom  aluminum 
arc  plates  with  some  of  the  CTI  ECAY  HR+  detector  modules  mounted  in 
place.  Each  detector  module  has  an  aluminum  guide  that  slides  along  a 
groove  in  the  arc  plate  and  is  stopped  (and  secured  against  the  arc  plate) 
using  a  mortise  and  tenon.  The  detector  modules  are  individually  angled  to 
point  towards  the  center  of  the  camera. 


Fig,  1:  ["Cjcholine  irmgc  of  prostate  cancer  (a)  before  and  (b)  after 
treatment.  These  grayscale  images  indicate  a  high  (white)  uptake  in  the 
prostate  cancer  compared  with  a  low  (gray)  uptake  elsewhere.  Images  are 
provided  by  Hara  and  co-workers  [8]. 

(a) 


possible  to  the  prostate,  which  improves  sensitivity  and  minimizes 
attenuation. 
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Fig.  5  3-D  drawing  of  the  prostate  PET  camera  with  some  gantr>'  details,  as 
viewed  from  the  side.  Only  a  few  detector  modules  are  shown  in  place. 
The  arc  plates  (with  detector  module  grooves,  shown  in  Fig,  4b)  are 
mounted  onto  thicker  aluminum  arc  back-plates  that  also  hold  the  lead 
shielding  in  place.  The  arc  plates  (white),  arc  back-plates  (dark  gray)  and 
lead  shielding  (light  gray)  are  symmetric  (-left  and  right  in  drawing)  since 
there  are  two  axial  rows  of  detector  modules.  The  upper  bank  is  mounted 
onto  a  main  aluminum  back-plate  in  order  to  control  the  vertical  motion 
using  a  hand-crank,  allowing  upward  movement  for  patient  access  and 
downward  to  maximize  sensitivity.  A  second  hand-crank  is  used  to  adjust 
the  overall  tilt  of  both  upper  and  lower  banks.  The  near  end  of  the  drawing 
shows  the  electronics.  The  six  DHI  boards  arc  mounted  horizontally  to  the 
housings  that  hold  their  corresponding  Analog  Subsection  boards.  The 
Coincidence  Processor  board  is  mounted  onto  the  cabinet  door.  The  power 
supplies  arc  mounted  on  the  far  end  of  the  gantry.  Most  of  the  cabling  is  not 
shown,  as  well  as  the  patient  bed. 


Fig.  6:  Statistical  noise-free  reconstructed  image  of  a  point  source  grid 
when  arcs  are  separated  by  the  ‘nominal’  distance.  The  point  sources  are  5 
cm  apart,  with  the  grid  covering  the  majority  of  the  field  of  view.  Most 
point  sources  are  located  between  pixels  {i.e.,  the  center  point  source 
contributes  to  four  nearby  pixels).  All  point  sources  are  visible  (black  = 
high  uptake).  Increased  radial  blurring  is  observed  near  the  outer  surface 
of  the  patient  as  expected  due  to  penetration  effects.  However,  good 
resolution  is  demonstrated  in  the  center  near  the  prostate,  with  minimal 
blurring  in  a  10  cm  diameter  central  region.  The  central  point  has  both  a 
radial  and  tangential  fwhm  of  2  mm.  The  surrounding  eight  points  have  a 


radial  fwhm  of  4  mm  and  tangential  fwhm  of  3  mm.  Improved  resolution  is 
expected  when  the  reconstruction  corrects  for  detector  penetration. 


Fig.  7:  Statistical  noise-free  reconstructed  image  of  a  point  source  grid 
when  arc  gap  distance  is  increased  by  +5  cm  from  ‘nominal’  position.  The 
point  sources  are  5  cm  apart,  with  the  grid  covering  the  majority  of  the  field 
of  view.  All  point  sources  are  clearly  visible  (black  =  high  uptake),  Good 
resolution  is  still  demonstrated  in  the  center  near  the  prostate,  with  minimal 
blurring  in  a  10  cm  diameter  central  region. 


(a) 


Fig.  8:  (a)  Body  outline  of  patient  with  2.5  cm  diameter  circular  “prostate” 
divided  into  quadrants  with  “tumors”  placed  in  the  upper  right  and  lower 
left  quadrants  (white  =  high  uptake),  assuming  a  tumor  to  background  ratio 
of  2; I.  The  “patient”  is  37.5  cm  wide  and  18.0  cm  deep,  (b)  Reconstructed 
image  of  the  phantom  using  an  iterative  maximum  likelihood  algorithm  that 
assumes  a  30%  scatter  fraction  and  20%  randoms  fraction.  Image 
represents  745  kcounts  in  a  single  slice,  which  should  be  achievable  with  a 
6  minute  scan  after  a  10  mCi  injection.  Both  tumors  are  clearly  visible. 
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Abstract — The  recent  development  of  new  prostate  tracers  has 
motivated  us  to  build  a  low  cost  PET  camera  optimized  to  image 
the  prostate.  Coincidence  imaging  of  positron  emitters  is  achieved 
using  a  pair  of  curved  detector  banks.  The  bottom  bank  is  fixed  be¬ 
low  the  patient  bed,  and  the  top  bank,  which  is  above  the  patient, 
moves  upward  for  patient  access  and  downward  for  maximum  sen¬ 
sitivity.  In  this  paper,  we  study  the  design  of  septa  for  the  prostate 
camera  using  Monte  Carlo  simulations.  The  system  performance 
is  measured  by  the  detectability  of  a  prostate  tumor  and  the  con¬ 
ventional  noise  equivalent  count  (NEC).  We  have  studied  27  septa 
configurations.  The  results  show  that  the  septa  design  has  a  large 
impact  on  the  lesion  detection  at  a  given  activity  concentration.  At 
the  background  activity  level  of  0.1  /xCi/cc,  sparse  septa  with  8- 
crystal  spacing  outperforms  the  traditional  two-dimension  (inter¬ 
plane  septa)  and  three-dimension  (septaless)  designs  in  terms  of 
both  the  lesion  detection  and  NEC.  Significant  differences  are  also 
observed  between  the  lesion  detectability  and  NEC  performance, 
indicating  that  the  NEC  does  not  correlate  well  with  this  lesion  de¬ 
tection  task. 


I.  Introduction 

Prostate  cancer  has  a  prevalence  and  diagnostic  rate  simi¬ 
lar  to  breast  cancer,  with  360,000  new  cases  diagnosed  each 
year  and  two  million  men  affected  by  the  disease  in  the  United 
States.  Prostate  cancer  suspicion  is  typically  based  on  an  ele¬ 
vated  prostate-specific  antigen  (PSA)  level  or  a  suspicious  node 
found  during  a  digital  rectal  exam.  Serum  PSA  values  do  not 
always  correlate  well  with  clinical  diagnosis  or  outcomes.  Pal¬ 
pation  is  subjective,  insensitive  and  inexact.  A  new  imaging 
technology  for  sensitive  detection  of  early  stage  prostate  cancer 
is  needed  to  confirm  initial  diagnosis  and  help  guide  treatment 
decisions. 

The  recent  development  of  new  prostate  tracers  has  shown 
promising  results  in  detecting  prostate  cancer  using  carbon-1 1- 
choline  [1].  This  has  motivated  us  to  build  a  low  cost  PET  cam¬ 
era  optimized  to  image  the  prostate  [2].  Coincidence  imaging 
of  positron  emitters  is  achieved  using  a  pair  of  curved  detector 
banks.  Fig.  1  shows  the  transaxial  view  of  the  camera.  The 
bottom  bank  is  fixed  below  the  patient  bed,  and  the  top  bank 
moves  upward  for  patient  access  and  downward  for  maximum 
sensitivity.  Each  bank  consists  of  two  axial  rows  of  20  block  de¬ 
tectors  (8x8  arrays  of  4.5  x  4.5  x 30mm^  BGO  crystals),  forming 
two  arcs  with  a  minor  axis  of  45  cm  and  major  axis  of  70  cm. 
Our  prostate  camera  has  about  one-fourth  the  number  of  detec¬ 
tors  as  in  a  conventional  PET  system  because:  (a)  the  patient 
is  not  fully  encircled  in  2D,  (b)  the  patient  port  is  smaller,  and 
(c)  the  axial  coverage  is  only  8  cm.  However,  since  the  average 
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distance  to  the  detectors  is  approximately  one-half  that  of  a  con¬ 
ventional  2D  PET  system,  the  solid  angle  for  a  central  source  is 
approximately  double  for  average-sized  men.  Thus,  we  expect 
to  achieve  improved  detection  efficiency  at  a  lower  cost. 

Individual  detector  modules  are  angled  to  point  towards  the 
camera  center  near  the  prostate  location  in  order  to  reduce  pen¬ 
etration  effects  (in  the  detector)  for  annihilation  photons  orig¬ 
inating  in  the  prostate.  This  maintains  high  resolution  selec¬ 
tively  in  the  region  of  the  prostate  and  is  an  unique  feature  of 
this  non-circular  camera  geometry.  Annihilation  photons  from 
other  parts  of  the  field  of  view  (FOV)  will  suffer  increased  pene¬ 
tration  effects,  but  these  regions  are  less  important.  Both  detec¬ 
tor  banks  can  be  tilted  to  image  the  prostate  while  minimizing 
attenuation  (i.e.,  above  the  buttocks  and  below  the  stomach),  but 
the  gantry  allows  zero  tilt  for  thin  patients.  Despite  this  special 
geometry,  we  are  able  to  reconstruct  nearly  artifact  free  images 
around  the  prostate  region  by  using  an  iterative  reconstruction 
algorithm  [2]. 

The  camera  design  includes  septa  to  reduce  background 
events  from  random  coincidences  and  from  Compton  scatter 
photons  in  the  patient.  Here  we  study  the  effect  of  different 
septa  designs  on  the  system  performance  using  Monte  Carlo 
simulations.  Septa  designs  have  been  studied  for  conventional 
PET  scanners  [3],  [4],  [5],  [6],  [7],  [8],  as  well  as  the  coinci¬ 
dence  gamma  cameras  [9],  [10],  [11].  In  this  work  we  focus 
on  the  specific  application  and  special  geometry  of  the  prostate 
camera.  To  better  assess  the  image  quality,  we  measure  the  de¬ 
tectability  of  a  prostate  lesion  using  an  ideal  linear  observer,  in 
addition  to  the  noise  equivalent  count  (NEC)  that  was  primarily 
used  in  previous  studies.  We  demonstrate  the  limitation  of  NEC 
for  evaluating  the  prostate  camera  design. 


Fig.  1.  A  positron  tomograph  for  prostate  imaging  [1]. 
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11.  Figures  OF  Merit 

A  very  popular  figure  of  merit  used  to  evaluate  PET  scanner 
design  is  the  noise  equivalent  count,  which  is  defined  as 


NEC  = 


T  S  kR 


(1) 


that  of  NEC  in  (1).  The  major  difference  is  that  the  SNR  focuses 
only  on  the  LORs  that  intersect  with  the  lesion  and  the  signal  in 
each  LOR  is  individually  weighted  by  its  variance. 

When  ri{y)  is  normally  distributed,  the  SNR  is  related  to  the 
area  under  the  ROC  curve  by  the  error  function  [13].  Therefore, 
we  use  (4)  to  measure  the  lesion  detectability. 


where  T,  5,  and  R  are  total  number  of  true  coincidences,  scat¬ 
tered  events,  and  random  events,  respectively.  Here  we  use 
k  =  1  assuming  the  randoms  are  not  pre-subtracted  and  the 
expectation  of  the  randoms  can  be  estimated  separately.  NEC 
represents  the  overall  signal-to-noise  ratio  (SNR)  of  the  data  set 
and  is  not  specific  to  any  local  region. 

Since  our  camera  is  designed  specifically  for  imaging  prostate 
tumors,  we  believe  that  it  is  more  appropriate  to  measure  sys¬ 
tem  performance  for  this  particular  detection  task.  We  choose 
a  “signal-known-exactly,  background-known-exactly”  (SKE- 
BKE)  task  because  the  location  of  the  prostate  is  relatively  fixed 
and  the  background  uptake  is  fairly  uniform  [1].  We  use  a 
prewhitening  (PW)  numerical  observer,  which  is  the  ideal  ob¬ 
server  under  Gaussian  noise,  to  evaluate  the  lesion  detectability. 
The  performance  of  the  PW  observer  measures  the  information 
content  in  the  data  for  lesion  detection. 

Let  hi  and  ho  denote  the  expectations  of  the  measurements 
with  and  without  a  lesion,  respectively.  The  PW  observer  com¬ 
putes  the  following  test  statistic  t)  for  a  given  data  set  y 

=  (2) 


where  S  is  the  ensemble  covariance  matrix  of  y.  A  decision 
whether  there  is  a  lesion  or  not  is  then  made  by  comparing  this 
statistic  to  a  pre-selected  threshold.  The  detectability  of  the  le¬ 
sion  can  be  measured  by  the  SNR  of  r?(2/),  which  is  defined  as 
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PET  data  can  be  modeled  as  independent  Poisson  random 
variables  with  the  expectation  y  related  to  the  tracer  distribution 
X  through  an  affine  transform 


y  =  Px  -f-  5  -fr, 


where  P  is  the  projection  probability  matrix  with  the  (i,  j)th 
element  denoting  the  probability  of  detecting  an  event  from  the 
jth  voxel  by  the  ith  LOR,  s  and  r  denote  the  expectation  of  the 
scattered  and  random  events,  respectively.  Using  this  model,  we 
can  show  that 
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where  f  is  the  tracer  uptake  of  a  prostate  lesion.  We  have  shown 
in  [12]  that  the  SNR  computed  in  (4)  is  the  same  as  the  SNR  of 
a  PW  observer  applied  to  images  reconstructed  using  the  max¬ 
imum  a  posteriori  (MAP)  principle,  but  not  using  the  filtered 
backprojection  method.  The  form  in  (4)  is  somewhat  similar  to 


III.  Monte  Carlo  Simulation 

We  use  Monte  Carlo  simulation  to  estimate  the  true  coinci¬ 
dence  rate,  scatters  and  randoms  under  each  septa  configura¬ 
tion.  The  Monte  Carlo  simulation  is  conducted  using  the  Sim- 
SET  software^ .  Both  Compton  and  photo-electric  interactions 
are  simulated  in  the  patient  body  and  septa.  The  patient  body  is 
modeled  as  an  elliptical  region  with  major  and  minor  axes  being 
40  cm  and  32  cm,  respectively,  which  is  the  medium  size  of  six 
randomly  selected  male  persons  in  our  department.  It  has  uni¬ 
form  activity  and  the  attenuation  coefficient  of  water.  The  axial 
length  of  the  body  in  the  simulation  is  60  cm,  beyond  which  the 
activity  has  little  contribution  to  the  data. 

The  current  version  of  SimSET  is  limited  to  circular  ring  PET 
scanners.  To  study  the  special  geometry  of  the  prostate  camera, 
we  model  the  septa  as  part  of  the  attenuation  object  as  shown 
in  Fig.  2.  The  material  of  the  septa  is  tungsten.  The  shape  of 
the  septa  is  modeled  as  the  difference  of  two  elliptical  regions 
with  major  and  minor  axes  of  the  inner  one  being  57  cm  and  32 
cm,  respectively,  and  those  of  the  outer  one  (1  -F  a)  57  cm  and 
(1  +  a)32  cm,  respectively,  where  a  is  a  constant  ranging  from 
0.04  to  0.64.  We  use  longer  septa  near  the  detectors  that  are  fur¬ 
ther  away  from  the  radial  center  to  limit  their  singles  rates.  We 
also  studied  two  septa  designs  with  uniform  length  and  found 
comparable  results.  The  top  and  bottom  axial  planes  contain 
perfect  absorber  to  prevent  photons  from  escaping  from  the  ax¬ 
ial  openings  so  that  the  acceptance  of  a  photon  is  completely 
determined  by  the  septa.  Therefore,  we  can  use  a  ring  detector 
surrounding  the  whole  object  and  the  resulting  sensitivity  is  the 
same  as  that  of  the  prostate  camera.  To  reduce  the  effect  of  the 
noise  in  Monte  Carlo  results,  we  rebin  the  data  into  one  transax- 
ial  sinogram  when  calculating  the  lesion  detectability  using  (4). 
The  thickness  of  the  septa  varies  from  0.4  mm  to  2.4  mm.  The 
axial  spacing  between  the  septa  varies  from  half  a  crystal  (2.4 
mm)  to  16  crystals  (77.6  mm).  Table  I  summaries  the  septa  de¬ 
signs  that  we  simulated. 

True  coincidence  events  and  scattered  events  are  binned  sep¬ 
arately.  For  all  the  results  presented  in  this  paper  we  use  an  en¬ 
ergy  threshold  of  355  keV,  except  where  noted  otherwise.  The 
singles  rates  are  estimated  by  running  the  simulation  in  SPECT 
mode  with  the  same  number  of  photons.  Random  events  are 
estimated  from  the  singles  rates  by 

rjt  =  r  •  singles(z)  •  singles(i),  (5) 

where  singles(z)  and  singles(j)  are  the  singles  rates  of  the  two 
detectors  forming  the  kih  LOR,  and  r  =  12ns  is  the  coincidence 
timing  window.  All  events  are  then  corrected  for  deadtime  as 
described  in  Section  IV. 

^Available  at  http://depts.washington.edur simset/html/simset_main.html 
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TABLE  I 

Parameters  of  the  simulated  septa  designs 


No. 

Spacing 

(crystal) 

thickness 

(mm) 

a 

1 

0.5 

MSM 

BH 

09 

BH 

0.08 

0.08 

5 

1 

0.4 

0.16 

6 

0.4 

0.32 

7 

0.8 

0.08 

8 

0.8 

0.16 

9 

0.8 

0.32 

10 

1.2 

0.08 

11 

1.2 

0.16 

12 

1.6 

0.08 

13 

1.6 

0.16 

14 

2.4 

0.04 

15 

2 

OH 

16 

17 

18 

4 

19 

1.2 

0.16 

20 

1.2 

0.32 

21 

1.6 

0.32 

22 

8 

1.6 

0.64 

23 

2.4 

0.16 

24 

2.4 

0.32 

25 

2.4 

0.64 

26 

16 

0.32 

27 

(no  septa) 

0.64 

IV.  Deadtime  MODEL 

The  deadtime  of  the  camera  is  modeled  similar  to  that  de¬ 
scribed  in  [14].  Two  parts  of  the  deadtime  are  considered,  one 
from  the  processing  of  single  events  in  each  detector  block  and 
the  other  from  the  coincidence  processing  circuit. 

For  each  detector,  the  live  time  can  be  modeled  as 

di  =  (6) 

where  N  is  the  number  of  detectors  per  detector  block, 
singles  (i)  is  the  average  singles  rate  of  the  individual  detectors 
within  the  block,  and  uiock  is  the  time  constant  of  the  detector 
block.  The  sensitivity  of  the  kth  LOR  due  to  the  detector  block 
deadtime  is  then 

Tik  =  di  ^  dj ,  (7) 

where  i  and  j  are  the  indices  of  the  two  detectors  forming  the 
kih  LOR. 

The  system  wide  coincidence  processing  live  time  is 
Coinc.  Live  = 

where  Tsys  is  the  time  constant,  and  Caine.  Load  is  the  ideal 
count  rate  experienced  by  the  coincidence  processing  circuit.  It 
is  calculated  as 

Coinc.Load  =  ^(4  +  5^:  -f  2rk)nk,  (9) 

k 


(b) 


Fig.  2.  Transaxial  (a)  and  coronal  (b)  views  through  the  center  of  the  attenuation 
image  used  in  the  Monte  Carlo  simulation.  The  inner  elliptical  region  is  the 
patient  body.  The  black  band  is  the  septa.  The  top  and  bottom  axial  planes 
contain  perfect  absoiber  to  prevent  photons  from  escaping  from  the  axial 
openings.  The  parameters  of  septa  shown  in  this  figure  are  2-crystal  spacing, 
1.2-mm  thick,  and  a  =  0.16. 

where  s^,  and  are  the  ideal  count  rates  of  the  trues,  scat¬ 
ters,  and  randoms  in  the  A:th  LOR,  respectively.  The  factor  of  2 
accounts  for  the  randoms  in  both  prompt  and  delayed  windows^. 
The  overall  sensitivity  of  the  kth  LOR  due  to  the  deadtime  loss 
is  Uk  *  Coinc.Live. 

V.  Normalization 

A  20  cm  diameter,  19  cm  long  uniform  cylinder  with  activ¬ 
ity  concentration  of  0.06  /iCi/cc  was  scanned  in  an  ECAT  HR 
scanner  (CTI  PET  systems,  Knoxville,  TN)  in  3D  mode.  The 
measured  singles  rates  are  compared  with  that  obtained  from 
the  Monte  Carlo  simulation  of  the  same  cylinder  scan  to  calcu¬ 
late  the  detector  efficiency.  After  scaling  the  simulation  result 
with  this  calculated  detector  efficiency  (=0.49),  the  true  coinci¬ 
dence,  scatters,  and  randoms  are  compared.  Fig.  3  shows  the 
measured  and  simulated  projection  profiles.  The  simulated  scat¬ 
ter  profile  matches  very  well  with  the  real  measurement.  The 
slight  discrepancy  between  the  measured  and  simulated  true  co¬ 
incidence  between  0  and  -100  mm  offset  is  due  to  the  attenuation 
of  the  patient  bed,  which  is  not  modeled  in  the  simulation.  The 
comparison  between  the  randoms  calculated  from  singles  rates 
and  that  measured  using  delayed  window  shows  a  nearly  perfect 
match  (a  ratio  of  1.0).  The  simulated  event  rate  of  trues  plus 
scatters  also  matches  the  measured  event  rate  well  (a  ratio  of 
0.97).  In  addition  to  the  detector  efficiency,  we  also  model  the 

^Here  we  assume  that  the  events  in  the  delayed  window  will  be  used  to  esti¬ 
mate  the  expectations  of  the  random  events. 


^  ft 


4 


missing  data  caused  by  the  gap  between  the  two  detector  banks 
(see  Fig.  1)  in  the  prostate  camera  simulations. 


Fig.  3.  Comparison  between  the  measured  ('+’)  and  the  simulated  (solid  line) 

projection  profiles. 

To  estimate  the  time  constants  rtiock  and  Tgys  in  the  dead¬ 
time  model,  another  uniform  cylinder  (the  same  geometry)  with 
activity  concentration  of  0.6  //Ci/cc  was  scanned  in  the  ECAT 
HR  scanner.  The  system  deadtime  as  reported  by  the  scanner  is 
about  23%.  The  fitted  time  constants  are  niock  -  2.2  fis  and 
'^sys  “  118  ns. 

VI.  Results 

A.  Lesion  detectability  V5.  NEC 

We  studied  27  different  septa  configurations  (Table  I).  Fig.  4 
shows  the  lesion  detectability  (SNR)  and  NEC  of  all  the  septa 
designs  at  different  background  activity  levels.  The  simulated 
lesion  is  a  2.5-cm  diameter  sphere  located  at  the  center  of  the 
FOV.  The  lesion  to  background  activity  ratio  is  2:1.  The  scan 
duration  is  1  minute.  Both  lesion  detectability  and  NEC  show 
that  longer  and  thicker  septa  are  preferred  for  high  activity  con¬ 
centration  due  to  the  reduction  of  scatters  and  randoms.  How¬ 
ever,  doubling  the  length  often  has  a  stronger  effect  than  dou¬ 
bling  the  thickness.  The  lesion  detectability  curves  clearly  show 
that  sparse  septa  (more  spacing)  can  improve  lesion  detection 
with  the  best  peak  SNR  achieved  by  design  no.  25  around  1.0- 
/iCi/cc  background  activity  level.  In  comparison,  using  sparse 
septa  only  increases  the  NEC  at  low  activity  levels,  but  not  the 
peak  NEC.  Other  significant  differences  are  also  observable  be¬ 
tween  the  SNR  plots  and  NEC  plots,  indicating  that  NEC  does 
not  correlate  well  with  this  lesion  detection  task. 

For  the  application  of  the  prostate  camera,  we  expect  that 
the  typical  background  activity  concentration  will  be  around  0.1 
juCi/cc,  which  corresponds  to  a  uniform  uptake  of  7.5-mCi  in¬ 
jection  into  a  75-kg  patient.  In  Fig.  5,  we  plot  the  SNR  and  NEC 
at  0.1  ^Ci/cc  for  all  the  septa  designs.  The  results  show  that  de¬ 
sign  no.  24  outperforms  others  in  terms  of  both  lesion  detection 
and  NEC,  and  so  is  the  best  septa  configuration  among  all  that 
we  have  studied. 


Design  No. 

(b) 


Fig.  5.  Performance  of  all  the  septa  designs  at  0.1/iCi/cc.  (a)  lesion  detectabil¬ 
ity;  (b)  NEC.  The  design  number  refers  to  Table  I.  The  dashed  lines  indicate 
the  change  in  septa  spacing. 

B,  Energy  resolution 

We  studied  the  effect  of  energy  resolution  on  septa  design 
by  varying  the  energy  threshold  in  the  Monte  Carlo  simulation. 
Fig.  6  shows  the  performance  of  a  system  with  perfect  energy 
resolution  where  we  set  the  energy  threshold  to  511  keV,  and 
Fig.  7  shows  the  performance  of  a  system  with  poor  energy  res¬ 
olution  where  the  energy  threshold  is  200  keV.  Comparing  the 
two  plots  at  the  same  activity  level,  we  found  that  good  energy 
resolution  increases  the  lesion  detectability  and  favors  less  septa 
due  to  the  better  rejection  of  events  scattered  inside  the  patient 
body  and  the  septa. 

VII.  Conclusion 

We  have  conducted  Monte  Carlo  simulations  to  study  the 
effect  of  septa  configuration  on  lesion  detection.  The  results 
clearly  show  that  properly  designed  sparse  septa  can  improve 
lesion  detection  over  the  traditional  two-dimension  (inter-plane 
septa)  or  three-dimension  (septaless)  configurations.  Among  all 
the  designs  that  we  studied  here,  we  found  that  the  septa  with 
8-crystal  spacing  give  the  best  lesion  detection  and  NEC  at  the 
background  activity  level  of  0. 1  jiCMcc.  In  future  work  we  will 
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(b) 


Fig.  4.  Performance  of  27  septa  designs  as  a  function  of  background  activity  levels,  (a)  lesion  detectability  plots;  (b)  NEC  plots.  The  legend  refers  to  the  design 
number  shown  in  Table  I. 


refine  the  search  for  the  optimal  septa  design.  We  will  also  con¬ 
struct  optimal  septa  designs  and  evaluate  their  performance  us¬ 
ing  real  measurements. 
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Fig.  6.  Performance  of  a  system  with  perfect  energy  resolution  (energy  thresh¬ 
old  =  51 1  keV).  (a)  lesion  detectability;  (b)  NEC.  The  design  number  refers 
to  Thble  I. 


MAP  reconstruction  using  computer  observers,”  IEEE  Trans  Med  Im,  vol. 
20,  pp.  815-822,  2001. 

[13]  H.  H.  Barrett,  G.  K,  Abbey,  and  E.  Clarkson,  “Objective  assessment  of 
image  quality.  III.  ROC  metrics,  ideal  observers,  and  likelihood-generating 
functions,”  Journal  of  Optical  Society  of  America  A,  vol.  15,  pp.  1520- 
1525,  1998. 

[14]  L.  Eriksson,  K.  Wienhard,  and  M.  Dahlbom,  “A  simple  data  loss  model  for 
positron  camera  systems,”  IEEE  Trans.  Nucl.  Sci,  vol.  41,  pp.  1566-1570, 
1994. 


(a) 


(b) 


Fig.  7.  Performance  of  a  system  with  poor  energy  resolution  (energy  threshold 
=  200  keV).  (a)  lesion  detectability;  (b)  NEC.  The  design  number  refers  to 
Table  I. 


